ZEROS OF WEAKLY HOLOMORPHIC MODULAR FORMS OF LEVELS 

2 AND 3 



SHARON ANNE GARTHWAITE AND PAUL JENKINS 

Abstract. Let M?(iV) be the space of weakly holomorphic modular forms for To(A^) that 
are holomorphic at all cusps except possibly at oo. We study a canonical basis for M|(2) 
and M^(3) and prove that almost all modular forms in this basis have the property that 
the majority of their zeros in a fundamental domain lie on a lower boundary arc of the 
fundamental domain. 



1. Introduction 

In studying a complex-valued function, it is natural to attempt to locate the zeros of the 
function; in fact, one of the most famous unsolved problems in mathematics asks whether 
the nontrivial zeros of the Riemann zeta function ((s) lie on a particular line. In this paper, 
we study the locations of the zeros of certain modular forms, and show that most of their 
zeros in a fundamental domain occur on a particular circular arc. 

For the Eisenstein series, perhaps the easiest examples of modular forms, a great deal is 
known about the locations of the zeros. In the 1960s, Wohlfahrt [22] showed that for even 
4 < k < 26, all zeros of the Eisenstein series E^{z) in the standard fundamental domain for 
SL 2 (Z) lie on the unit circle \z\ = 1, and R. A. Rankin [T7] extended the range of values of k 
for which this holds. Shortly afterward, F.K.C. Rankin and Swinnerton-Dyer [16] proved this 
result for all weights k > 4. R.A. Rankin [18j also obtained the result for certain Poincare 
series, which generalize Eisenstein series. Similar results have been obtained for Eisenstein 
series for Tq(2) and Tq(3) by Miezaki, Nozaki, and Shigezumi [H], for Eisenstein series for 
To(5) and 1^(7) and for Poincare series for rj$(2) and Tq(3) by Shigezumi [201 121] , and for a 
family of Eisenstein series for T(2) by the first author, Long, Swisher, and Treneer [XT] . 

The above results which locate the zeros of Eisenstein series and Poincare series use the 
same general idea of approximating the modular form by an elementary function having the 
required number of zeros on the arc. For example, the Eisenstein series Ej-(z) may be written 
as a sum over an integer lattice. When z is restricted to the unit circle, so that z = e 10 , the 
four main terms of this series combine to give 2e~^~ cos (^) . Rankin and Swinnerton-Dyer's 
proof shows that the additional terms are small, so the zeros of Ek{z) are close to the zeros 
of this trigonometric function. 

In 1997, Asai, Kaneko, and Ninomiya [3] used this idea to study the zeros of polynomials 
related to the modular function j(r). The j{r) function is a Hauptmodul, or an isomorphism 
from the quotient of the upper half plane EI under the action of SL 2 (Z) to the complex plane 
C. It generates all modular functions on SL 2 (Z), and it also parameterizes the isomorphism 
classes of elliptic curves over C. The image of j(z) under the Hecke operator T n is a polyno- 
mial in j(z), which we write as P n (j(z)). Letting A(z) be the modular discriminant, which 
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is a weight 12 cusp form on SL 2 (z) with no zeros in H, a generating function for the P n (j(z)) 
is given by 

Using this generating function, Asai et al. approximated the polynomials P n (j(z)) by trigono- 
metric functions well enough to prove that their zeros in the fundamental domain lie on the 
unit circle. 

Duke and the second author [5] extended the results on j(z) to a two-parameter family of 
modular forms that form bases for spaces of weakly holomorphic modular forms of level 1. 
In this case, the connection to elementary functions is less direct. Cauchy's integral formula 
relates the modular forms to a contour integral of a generalized version of the generating 
function (1.1). An application of the residue theorem produces the elementary functions, 
and the proof follows by bounding the integral over a range of values for r and z. The zeros 
again lie on the unit circle for many of the forms in the family, though in contrast to previous 
results, it is known that this property does not hold for all of the modular forms. 

We mention one further result using a different technique by Hahn [T3], who obtained 
general results on the zeros of Eisenstein series for genus zero Fuchsian groups; the general 
idea is an analogue of the classical argument that shows that the zeros of an orthogonal 
polynomial all lie on the real line. 

This question of locating zeros of modular forms is made even more interesting by results 
of Rudnick [19] that showed that the zeros of Hecke eigenforms of weight k, in a sense the 
orthogonal complement of the Eisenstein series, are expected to become equidistributed in 
the fundamental domain as k — > oo. Ghosh and Sarnak [12] gave a lower bound for the 
density of zeros lying on certain arcs for such eigenforms. In a different direction, Basraoui 
and Sebbar [B] proved that the quasi- modular form E 2 (t) has infinitely many zeros that are 
inequivalent under SL^Z), and that none of these lie within the fundamental domain. 

In this paper, we examine modular forms in a basis for certain spaces of weakly holomor- 
phic modular forms of arbitrary integral weight and levels 2 and 3. We show that for almost 
all of the basis elements, most of their zeros in a fundamental domain for To (2) or To (3) 
lie on a circular arc along the lower boundary of the fundamental domain. This is possible 
because we can again approximate these modular forms by elementary functions; however, 
the shape of the fundamental domain makes it difficult to accurately locate all of the zeros. 



2. Definitions and statement of results 

Let Mfc(2) be the space of holomorphic modular forms of weight k for the group To (2) = 
{(cd) e SL 2 (Z) : c = (mod 2)}, and let M[{2) be the corresponding space of weakly 
holomorphic modular forms, or modular forms that are holomorphic on the upper half plane 
and meromorphic at the cusps. Let M|(2) be the subspace of M\(2) consisting of forms which 
are holomorphic away from the cusp at oo. This space appears, for instance, in [15], where it 
is shown that traces of negative integral weight forms in such a space appear as coefficients 
of certain half integral weight forms of level AN. Modular forms in M|(2) have been studied 
by Ahlgren PQ, who gave explicit formulas for the action of the ^-operator on forms in these 
spaces and obtained formulas for the exponents of their infinite product expansions, and by 
Andersen and the second author [2], who gave congruences for the coefficients of a basis for 
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For the group r (2), we use a fundamental domain in the upper half plane bounded by the 
lines Re(z) = — \ and Re(z) = | and the circles of radius | centered at z — — | and z —\. 
We include the boundary on the left half of this fundamental domain, which is equivalent to 
the opposite boundary under the action of the matrices (J }) and (jj> ?)■ The cusps of this 
fundamental domain may be taken to be at 00 and at 0. 

2n 



1.5- 




Modular forms in M\{2) with real coefficients demonstrate a nice property on the lower 
boundary of this fundamental domain. For a Fourier series f(z) = a{n)e 2mnz with real 
Fourier coefficients a(n), note that f(a + bi) = a(n)e~ 2mna e~ 2nnb = f(—a + bi), or f(z) = 
f(—z). Thus, for modular forms / of weight k on r (2), if we let z = —\ + \e- 16 for < 6 < |, 
so that z is on the lower boundary of the (symmetric) fundamental domain, we find that 
(io) z= _L_ = i_ i e -i0 = _z. Therefore, f(z) = (2z + iy k f (^) = e~ ike J(z), and the 

normalized modular form e^r/(— \ + \e ie ) is real valued for 9 between and |. 

It is useful to define three particular modular forms of level 2. As usual, let q = e 2mz . Let 

m = (^y) 24 = V 1 - 24 + 276g + . . . G M« 

be the Hauptmodul for r (2). This form has integer coefficients, has a pole at 00, and 
vanishes at 0. Moreover, by the above argument <f>(z) is real-valued on the lower boundary 
of the fundamental domain, taking on values in [—64, 0]. The special value (— | + |) = —64 
arises from the relationships between <f>(z) and j(z) given by 

(f)(z) 2 4>(z) 

By using the fact that j (— \ + |) = = 1728, it is easy to see that (p {—\ + |) must be 
a common root of the polynomials (x + 256) 3 — 1728a; 2 = (x + 64) (x — 512) 2 and (x + 16) 3 — 
1728x = (x- 8) 2 (x + 64). 
Next, let 

F 2 (z) = 2E 2 (2z) - E 2 (z) = 1 + 24 ^ j ^ d g n 

n \d|n,dodd / 
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be the unique normalized holomorphic modular form of weight 2 and level 2. Here E 2 (z) 
is the weight 2 Eisenstein series E 2 (z) = 1 — 24 Yl^Li o~(n)q n . The form F 2 (z) has integer 
coefficients and a single zero at the elliptic point —\ + \- This can be seen by noting that 
i^zz+l) = F 2 {z + 1) = F 2 {z) at z = —\ + |, so an application of the modular equation 
yields —F 2 (z) = F 2 (z). Uniqueness comes from applying the valence formula for To (2) found 
in, for instance, [7]. Additionally, we define the Eisenstein series S^z) 6 M^{2) as 

(2-1) S*(z) = E ^~Z^ = q + 8q * + 28g 3 + ■ ■ ■ ; 

it is easily checked that S4 has integral Fourier coefficients and vanishes at 00. It does not 
vanish at the cusp at 0, as the valence formula shows that there are no cusp forms of weight 
4 and level 2. 

We now use these forms to construct a basis for M|(2). For general even weight k, we 
write k = 4£ + k', where k' G {0, 2}. A basis for M|(2) is given by 

for all integers n > —£. We note that this is an extension of the basis given in [4] for M^(2); 
similar sequences of modular forms for many levels appear in [S]. The basis elements are 
constructed as follows. We first define F — 1 and set fj?_i(z) = Sl(z)F k /(z). Next, for 

each n > — £ we define fj^l(z) inductively by multiplying /f^-iiz) by 4>{z) and subtracting 
off earlier basis elements. Note that since 5*4, Fy, and 4> have integral Fourier coefficients, 
each fkn( z ) has integral Fourier coefficients and is of the form Sl(z)F k >(z)F((f)(z)), where 
F(x) is a polynomial with integer coefficients of degree n + £ = n + |_f J ■ (In fact, F(x) is a 

generalized Faber polynomial; see [9], [ID].) Thus, if all of the zeros of the basis elements n 
lie on the lower boundary of the fundamental domain, then all of the zeros of the polynomial 
F{x) must lie in the interval [—64,0]. 

The main result of this paper is the following theorem. 

Theorem 1. Let f£l(z) be as above. If £ >0andn> U£ + 8, or if £ < and n > 15|£|+8 ; 

then at least [^n + |J of the n + |_|J nontrivial zeros of fjf^(z) in the fundamental domain 
for T (2) lie on the lower boundary of the fundamental domain. 

We note that the bounds n > 14£ + 8 and n > 15\£\ + 8 are not sharp, and that often 
many more of the zeros lie on the arc. In fact, for certain weights k close to 0, all of the zeros 
of all of the fj® (z) are on the lower boundary of the fundamental domain. However, some 

restriction on n in relation to £ is necessary, as there are also examples of f k 2 l(z) with zeros 
elsewhere. We discuss this further in Section [6} Additionally, we obtain similar results for 
a family of modular forms in M|(2) whose coefficients are dual to the Fourier coefficients of 

(2") tt 

fk n\ z ) an d f° r a basis for the space M^(3), showing that many of the zeros of these modular 
forms in the appropriate fundamental domain lie on the lower boundary. 

The remainder of this paper proceeds as follows: in Section [3] we give a generating function 
for the basis elements f k ^{z) and approximate their values on the lower boundary of the 
fundamental domain by a trigonometric function. In Section [4] we bound the error term to 
show that most of the zeros lie on the appropriate arc. Section [5] gives technical details on 
bounds for the error, and in Section [6] we discuss extensions of the main theorem to other 
modular forms for To (2) and to forms for To (3). 
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3. Generating functions and integration 

In this section we use Cauchy's integral formula to relate the basis elements f k n( z ) to a 
trigonometric function. Letting r = e 27nr , a generating function for the basis elements /^(z) 
is given in El-Guindy's paper (HI Theorem 1.2] as 

{ } hi "(^,)(r)0(r)-0(,)- 

Multiply by r~ n ~ x and integrate around r = 0. Changing variables from r to r and 
rewriting 0(t)F 2 (t) as a derivative, we obtain 



-i,2MA{Si4>F k '){r) 0(r)-0(z) -2ttz 



where A is some real number larger than |. We move the contour downward, noting that 
there is a pole whenever 0(r) — is zero, which happens when r is equivalent to z under 
the action of r (2). There are no other poles, since S4 has no zeros in the upper half plane, 
and if k' = 2, the zero of F 2 (r) is canceled by the zero of 4-(j>{r) = 0(r)F 2 (r). The closed 
contour that is the difference between the old integral and the new integral is in the clockwise 
direction, so we get a factor of —1 in Cauchy's integral theorem, and at each pole, we obtain 
a term of — 27fi times the residue of the integrand. 

If a function /(r) has a zero of order n at To, its logarithmic derivative has a simple pole 
with residue n. In calculating the residue, note that part of the integrand is the logarithmic 
derivative of 0(r) — <f)(z), which has a simple zero exactly at the values we are looking at, 
since is a Hauptmodul for To(2). This means that the logarithmic derivative of 0(r) — 4>(z) 
at points equivalent to z under r (2) will just give us a factor of 1 in the residue. Supposing 
that t = 72; = for some 7 G r (2), this is then multiplied by 



-2Trinyz 



{S{F k ,){z) 



-2m (SiF k ,)(jz) 



Since the denominator is a modular form of weight k on r (2), we have (S^F^^z) = 
(cz + d) k (SlF k i)(z) for 7 = (° b d ), and the residue becomes zj^i{ cz + d)~ k e~ 2nin ^ z \ 

Assume that z is on the lower boundary of the fundamental domain for r (2). The first 
two points where r is equivalent to z through which the contour moves are r = z and 
r = i^i- Calculating the residues, we find that 



fW (z) _ e -™» _ {2z + I)" V 2 ™^) - f MF„)(z) Ufa) - <K*» e-™ nT dT 
hA ) { + ) ~ Jc{S{F k ,){r) <Hr)-m -2m ^ 

where C is a contour that moves from left to right across the fundamental domain and passes 
below the points r = z and r = and above all other points equivalent to z under the 
action of r (2). 

We write z = — \ + \e ld for some < 6 < |, so that = \ — \ e ~ l °\ then the quantity 
e" 2 ™ 2 + {2z + l)- fc e " 2 ™(^i) can be simplified to 



6 SHARON ANNE GARTHWAITE AND PAUL JENKINS 

Putting all of this together and multiplying through by e lke /' 2 e -^ nsme ^ we enc [ U p with 

e "/V™-»7S [~\ + ~ (-l)"2cos - Trncos^ = 



ikd/2 — nn sin 8 2-iririT j_ 

' ' 7 c (^)(r)0(r)-^) e 

By the argument in Section [2j the left hand side is a real- valued function of 9. We note that 
the cosine function takes on alternating values of ±2 whenever ^ — nn cos 9 is equal to ran 
for m G Z. Since this quantity moves from — ri7r at 6* = to ^ at 9 = ir/2, we know that 
there must be at least n + 1 + |_|J times that this happens. Thus, if we bound the integral 
term in absolute value by 2, then by the Intermediate Value Theorem we must have at least 
n + [|J zeros of the modular form on this arc. 

The dimension of the space of holomorphic modular forms M*.(2) is +1, and we get 
at most |_f J zeros not at elliptic points for these forms by the general valence formula in [7]. 
(Note that F 2 (z) has a zero at —5 + 5, which is already on the arc in question, so fj?l has a 
trivial zero there if k' — 2.) Counting a pole of order n at 00 and no other poles gives us a 
total of n+ |_f J zeros of the basis element fj^(z) G M\{2) whose locations are unknown. This 
argument proves that if the weighted modular form is close enough to the cosine function, 
then all of these zeros must be simple and must be on this arc on the lower boundary of the 
fundamental domain. 

Unfortunately, it is difficult to move the contour down far enough to prove that all of the 
zeros are on this arc; as the weight k or the order n of the pole increases, the contour will 
need to get closer and closer to r = if 9 is close to 0. Recall that for every fixed value of 
z, we need to prove that the integral is bounded by 2, after moving the contour below that 
value of z. As z gets close to the real line, this becomes very difficult-either the contour is 
not straight and the integral is harder to estimate, as r does not have a fixed imaginary part, 
or the contour must pass through more residues, adding additional terms to the equation. 

We can still prove that the majority of the zeros do indeed lie on this arc by choosing 
a fixed height for the contour, estimating the value of the integral along that contour, and 
showing that its absolute value is bounded above by 2. The goal then is to choose a contour 
low enough to capture as many zeros as possible, yet high enough to avoid additional residues 
and to avoid large values inside the integral. We choose r = u + | for \u\ < |, so that the 
contour has constant imaginary part |. 

For this choice to work, we must also limit the range of z, so that our contour passes below 
t = z and r = but above other images of z under r (2). If z = — \ + \e ie for | < 9 < |, 
then z has imaginary part > |, and a contour at a height of | picks up residues at r = z 
and t = ^fj-j- but no other points equivalent to z under the action of To (2); the maximum 
possible imaginary part of such a point is |. In this case, the quantity (4p — cos inside 
the cosine function has the value || — 7rn^ at 9 — | and the value at 9 = |, and passes 
through at least |_§ + multiples of n. Bounding the integral by 2 for the appropriate 
fkn( z ) wm finish the proof of Theorem 1 
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4. Bounding the integral 



In this section we bound 



e ike/2 e - 



(SjF k ,)(z) 0(r)F 2 (r) 
c(SiF k ,)(r) 4>(t) - <t>(z)' 



-2ninT 



dr 



for the values z = 



We will also 



\ + \e ie with 9 G [f , f] and r = u + \ with u G [-§, |] 
give some indication of how this bound might change if we allow 9 to approach and alter 
the countour accordingly. Details for the computation of the numerical bounds that appear 
here are provided in the next section. 



We seek a bound for 



e 2 e 



ik6 

e 2 e" 



i0f(2) 
J k,n 



2 2 



2 T 2 C 



A|0 



(-l) n 2cos 
^(«+i)i? 2 («+i) 



nn cos 9 ) = 



which is real- valued, by something less than 2. In absolute value, this integral is 



du, 



3 — 7rn(sin6— |) 



(52^)(-| + 



-2-ninu 



(SiF k ,)(u + 



<t>{u+l)-<t>{-- 2 + 
We have chosen Qffr 



du 



Consider the exponential term e 7rn ( sin6 ' 5). We have chosen ^(t) = |, so that sin 9 — | > 
for G [|, |], and this term has exponential decay as n — > 00; in this case e-^C 811161 -!) < 
.73041. If we find an upper bound for the absolute value of the integral, then for large enough 
n the right hand side is indeed less than 2, and we can apply the Intermediate Value Theorem 
as desired. It turns out that we can find a bound for the absolute value that removes the 

dependence on n, but may be exponential in t. However, if n is large enough in relation to 

(2) 

£, then we will see that for a fixed weight k, all but finitely many of the f k ni z ) have zeros 
on the appropriate arc. 

We note that the absolute value of the integral is certainly bounded above by 



Sa(- 



i + W0)F 2 (u+l)) 



<t>{u+i) 



2 T 2 C 



i.0 



du, 



F k >(u+l) 

and already the dependence on n has vanished. For the terms involving S4 and F k >, we find 
an upper bound for the maximum possible value of these terms over the appropriate ranges 
of u and 9, and pull these upper bounds outside of the integral. This leaves us with the 
contribution from 

0(u+i) 



<p(u 



du, 



which we consider in pieces. 
Consider first the quantity 



s 4 (- 



2 ~ 2 C 



SAu + 



Computations, explained in more detail in the next section, yield 



.014 < 



Sa 



u 



< 2.44141. 
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As a power series in q, we know that 5*4 has positive coefficients, as seen in (2.1) 
maximum value of \SAt 



The 



occurs when u — 0, when q is real and positive. Heuristically, the 
minimum value should occur when there is maximum cancelation between terms, or when 
q is real and negative, so that u = ±|, and we confirm this computationally. Decreasing 
S(r) both increases the upper bound, as we are adding larger positive terms, and potentially 
decreases the lower bound due to cancelation. 

Heuristically, the maximum value of |S4(z)| should occur when (—1 + cos 9) is close to 0, 
meaning 9 is close to 0, as here q is real, positive, and close to 1. Similarly, the minimum 
value of Si(z) should occur when (—1 + cos9) is close to ±1, meaning 9 is close to |; here 
q is real but negative, so there is extensive cancelation when adding terms. In this case, 



though, the size of q depends on 9, as we have |e 
the minimum indeed occurs at 9 



2m{ 



2 ■ 



where |e" 



-7r sin( : j) I 



2m[ 



-7r sin(#) I 



1 )\ = | e -™n(0)|_ Computationally, 
.04322. In general, we have 

< e- w/2 « .20788. 



We compute that 



.03 < 



Si 



- + -e l 
2 2 



< .99995. 



Moving the lower bound on 9 closer to increases the maximum value of |S , 4(^)|, though it 
does not appear to affect the minimum value. 
Putting this together, we have, for i > 0, 



SAu + 



and for i < 0, 



Next, we consider the term 



SAu + 



< 171.4251 



< [81.380341^1. 



Fk< 



1 _l ± p i0* 
2^2^- 



F 2 (u + 



F k ,(u + 



If k' = 2, this is \Fo 



which is bounded above by 8.00067. If k' = 0, this is 



\F 2 (u + |) , which is bounded above by 12.50005. Either way, the contribution is no more 
than 12.50005. Note that Fi(t) has positive coefficients, and is therefore large when SAt) 
is large. 



Finally, we consider 



¥ 6 ) 



du. 



The Hauptmodul <fi(z) is an injective mapping on the fundamental domain; it is real-valued 
and strictly decreasing on 9 G [0, |]. We restrict ourselves to 9 6 [|, |]. Computation shows 
that <j)(z) G [—64, —.033] on this domain. In contrast, <p(u + |) takes on a wide range of 
values, including some with very large and some with very small modulus. Bounding the 
numerator and denominator separately yields a trivial upper bound of roughly 79000, while 
numerical calculations indicate that the actual maximum is a little larger than 1. 
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In order to achieve a sharper bound, we will instead consider a related quantity, 



(4.1) 



where we use 6(z) to indicate 



D(z, 



T 



T 



2 ~ 2 C 



and 0(r) to indicate 4>{u + |) for ease of notation. 



This quantity is related to our Hauptmodul expression by the identity 



(4.2) 



i 



\1 + D(z,t) 



It is easier to work with D(z, r), as we know the numerator is real valued and within a small 
range, and bounding this quantity proves useful in Section [6] when discussing extensions of 
Theorem [TJ 

We will break our path of integration into pieces, and consider <J)(t) in relation to <fi(z) 
on each. As we know that <fi(z) is real, we consider the real and imaginary parts of 0(r) = 
<p [u + |) separately for u G [— |, |]. It is clear that e 2mT = (cos(27ru) + i sin(27rw)), and 
so 9?(0(u + §)) = + §)), while 9f(0(« + §)) = -3f(0(-« + §)). With this in mind, 

we restrict our calculations to u E [— |, 0] and use symmetry for it e [0, |]. 

The numerical techniques described in the next section reveal that on the interval u G 
[-.5, -.21516], we have either 3£(0(r)) > 0, 3?(0(r)) < -128, or |3(0(r))| > 64; it follows 
that \D(z,t)\ < 1 on this interval. 

Next, we note that since cj){z) is real, then if we have the bound 3f(0(r)) > A, it follows 
that 

0(z) 



D(z,t)\ < 



(K(0(r))-0(z)) + Az 

If 3?(0(r)) > 0, then this is bounded above by 1, while if 9fJ(0(r)) < 0, then the maximum 
possible value of the right hand side for a fixed t as z varies is 



+ 1, 



2 +?R(^(r)) s 



. A lower bound for 5f(0(r)) thus gives us an upper bound 



occurring when <f>{z) - 
for \D(z,t)\. 

For u e [-.21516, -.18884] we have 3(0(r)) > 1 and K(0(r)) > -.0175. If 3?(0(r)) > 0, 
then our bound is 1, and if > 9ft(0(r)) > —.0175 then our bound is 1.00016. Either way, 
for u G [-.21516, -.18884], we have \D(z,t)\ < 1.00016. Similarly, on [-.18884, -.12878] 
we have S(0(r)) > .033, and we obtain a bound of 1.13192. 

Finally, we consider [-.12878, -0]. We have 9fc(0(r)) > -.01424, and so 



\D(z,t)\ < 



-.01424 



The maximum value occurs when </>(z) = —.033, and we find \D(z,t)\ < 1.75344 here. 

Altogether, we have \D(z, r)| < 1.75344. By breaking the integral into pieces, we compute 
more precisely that 

•1/2 

\D(z,r)\ du < 2 • 0.60496 = 1.20992. 

'-1/2 
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The relationship in (4.2) allows us to conclude that 



(u + 



i + ¥ 6 ) 



du < 2.20992. 



We can improve this further by noting that the fact that 9?(</>(t)) > on [—.45787, —.22531] 
implies that the integrand is bounded by 1 here, yielding 



u 



-du 



< 1.74520. 



We see that if we extend 9 closer to 0, and hence also decrease 3(r), this term has the 
most potential to blow up near u = 0, as this is where (p(z) and 9ft(0(r)) are both small. 
Additionally, if < 6 < |, then any contour with fixed imaginary part less than |, such as 
t = u + |, will cross the arc z = — \ + \e ld and at least one of its images under r (2), so 
restricting our 9 values is necessary to avoid a zero in the denominator. 

Putting all of these pieces together and using the fact that sin# is decreasing on [|, f], we 
see that for £ > 0, 



-7rn(sin6— |) 



{S{F k ,){u+\) </>(u+' 



-2ninu 



du 



< .73041 n |71.425| £ (12.50005)(l. 74520). 



Note that (.73041™) (12.50005) (1.74520) < 2 if n > 8, and (.73041 n )71.425 < 1 if n > 14; 
hence, the integral is less than our desired bound 2 if £ > and n > 14£ + 8. Similarly, 
for £ < 0, we replace |71.425| £ with |81.38034|' £ ', and find that our integral is bounded by 
2 if n > 15\£\ + 8. We can then apply the Intermediate Value Theorem to prove that the 
appropriate number of zeros are on the desired arc. 



5. Rigorously computing upper and lower bounds 

In the previous section, while bounding our integral we used upper and lower bounds on 
Eisenstein series and the Hauptmodul for values on a circular arc on the boundary of the 
fundamental domain and on a straight line segment. In this section we justify those bounds. 

It is useful for most of these calculations to truncate each series. For a modular form / 
with Fourier series / = Y2 a f( n ) < l n ^ we wm choose a positive integer iV and let / be the 
truncation of the Fourier series of / up to and including the q N term, and we let Rf = f — / 
be the remaining tail of the series. We bound / and Rf separately. 

The calculations for the Eisenstein series are straightforward, as we have explicit formulas 
for the Fourier coefficients in terms of divisor functions, while calculations for <fi require a 
little more finesse. We do not have a nice formula or a sharp bound for the growth rate 
of the Hauptmodul coefficients, and they are quite large, so there are more terms making a 
significant contribution to the value of the series. We begin by bounding the values of the 
Eisenstein series, and then use those bounds to tame d(z). 



By (2.1 ) we have 

oo 

S A (z) = J2 (*s(n) - *s (|))? n 
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For k > 1, we can generously bound <7fc(n) = ^2 d i n d 

< n h ¥+n k+1 . If 



pairs of divisors (<i, ^j. Thus, <7fc(n) 
by 



by yfn ■ y/n + y/n ■ n k by considering 
< t, then we can bound RS^z) 



-,2-Kiz I 



oo 

E 

n=N+l 



(crs(n)-a 3 (|)) 



< 



n=A r +10 

E 

n=N+l 

oo 

E 

n=AT+ll 



(n 2 + n 4 )t"- E (l + n 3 )t 2n . 

The last term reflects a lower bound on the summand 03 (|) coming from the even terms. 
Standard Taylor series methods involving derivatives of the geometric series (1 — x)^ 1 = x n 
taken at x — t allow us to bound each infinite series. For example, 

N+10 



E (» 2 

n=N+ll 



+ n 4 )t n 



t 2 + t 



t 4 + lit 3 + lit 2 + t 

f — ^-^ E 

n=l 



(i-ty 



n 



2wiz I 



< t 



-27rsin(f) 



-2tt 

e * . 



We find that 



We take N = 50. For S^z), we note that |e' 
I^WI < 2.86404- 1(T 23 . 

To find an upper bound for |<S4(;z)|, we explicitly compute X)n=i (°"3( n ) ~~ °"3 (f)) ^™ w ith 
t = e^r and then add the tail. We find that ^(-z)! < .99995 on the appropriate arc. 
For a lower bound on ^(z)!, we consider the real and imaginary parts of S^(— \ + \e l6 ) 
separately. We compute the first derivative with respect to 9 for each part, and for 9 e [|, |], 

each derivative has a trivial upper bound of Y^n=i 7Tn ' { a s( n ) ~ °3 (f )) (e _27r / 4 ) n < 8.01. A 
computation reveals that either the real or the imaginary part of S±{z) is at least .0302 when 



\e l ° with 9 



for all < n < 40000. The bounds on derivatives 



2 1 2 V6 1 3 40000' 

and on the tail RS^z) limit how close ^(z)! can get to 0, and we therefore conclude that 
IS^z)! is bounded below by .03. 

We can do similar upper bound calculations for S^r), F 2 (t), and F 2 (z), using the addi- 
tional fact that \e 2niT \ < t — e _27r/5 . For F 2 , we bound the tail by 



\RF 2 (t) 



E 

n=N+l 



n=7V+10 00 

ajinur < E a F2 (n)t n + 24 E Q ' (' 

n=N+l n=7V+ll 

where the use of | instead of n comes from the fact that we only consider odd divisors. We 
again take N = 50 in each case, and find that |S 4 (t)| < 2.44141, \F 2 (z)\ < 8.00067, and 
|F 2 (r)| < 12.50005. 

To compute a lower bound for [^(r)], we trivially bound the derivatives of the real and 
imaginary parts of S±{u+ |) by 

50 

E 2™ • a S4 (n) (e- 2n/5 ) n < 48.83. 

n=l 

We then compute the real and imaginary parts of S^u + |) for u = — \ + \ ■ 12 3 000 , where 
< n < 123000, verifying that at least one of these values is larger than .014010. We 
conclude that (^(r)] is bounded below by .014 on u e [— |, 0], and use symmetry to extend 
this to [0,^]. 
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We now consider the Hauptmodul 4>{j). In the previous section we needed information 
about the size of the real and imaginary parts of 0(r) and the value of <p(z). For these 
computations, we work with the truncations <p(z) and 0(t), taking into account the growth 
of the real and imaginary parts of the truncations and the error caused by ignoring the 
tail. In this case, the trivial bound on the partial derivatives is much larger than in the 
Eisenstein series case, so we truncate each series up to and including the g 30 term to shorten 
our computation time. 

We can express in terms of Eisenstein series of level 2 as 



4(z) 



E A {2z) 
S A (z) 



-16; 



we use this representation to bound R<f>(z). Observe that if we truncate (f)(z), then the tail 
satisfies 

R<P(z) = <f>(z) - fa) 
E A (2z 



Sa(z) 



- 16-0(2) 



(z) + 16 + ^2z)-(fa) + l6)S^z)) _ 16 _ , {z) 



S 4 {z) 



S 4 (z) 



< 



fa) + 16 



SAz) 



S A {z) - RS A {z) 



+ 



E A (2z) - (fa) + lQ)S A (z) 


+ \RE A {2z)\ 




S 4 {z) - RS A {z) 





We can now bound R<p(z). We compute bounds for the g 31 to q A0 terms directly to find that 
they contribute at most 6.46551- 10 -8 , and use the truncation formula to bound the remainder 
of the tail. In applying the formula, we truncate all Eisenstein series at N = 50 to use the 
previously computed upper and lower bounds. Note that, working as before, we can use 
the fact that E A (2z) = 1 + 240 ^f =1 (r 3 (n)q 2n to prove that \RE A (2z)\ < 6.40309 • 10~ 29 and 
\RE A (2r)\ < 2.16794- 10" 22 . We compute bounds for 4>(z) + 16 and E A (2z) - (4>(z) + 16)S A {z) 
by summing ^2\df{n)\e^^~ for each series. We find that the bounds are 544.01429 and 
7.29909 • 10~ 13 , respectively. Putting all of this together we have, for iV = 30, 



\R(f)(z)\ < 6.46551 • 10~ 8 + 544.01429 



.03 



7.29909 • 10" 



.03 

13 + 6.40309 • 10" 29 



7.05863 • 10" 30 



-1 + 



.03 - 7.05863 
< 6.46754 • 10" 8 . 



10 



-30 



We repeat the same calculations for 0(r). We find that 



\R(f>(r)\ < .001371 + 2593.07795 



.014 



.014- 2.86404- 10- 23 
6.40510 • 10~ 7 + 2.16794 ■ 10" 22 



1 + 



.014 - 2.86404 • 10" 23 
< 0.00142. 
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Now that we have bounds on the error caused by truncation, we numerically compute that 
+ \e ie ) < -.033 for 6 G [|, |] by calculating $(z) for iV = 30 at 9 = § and adding the 
bound for the tail \R(f>(z)\. 

The bounds for 4>(r) are slightly more difficult. As with lower bounds for the Eisenstein 
series, we consider the real and imaginary parts of 4>{r) separately. The maximum possible 
growth rate for the real part is 



d 

du 



30 



a^(n)e 5 cos(27rmt) 



< 



30 

2vm • \a^{n)\ (e^)" < 101197.78 



and the same bound holds for the derivative of the imaginary part. We again compute values 
of the real and imaginary parts of <f>(r) at a sampling of points, and use these bounds on the 
derivatives to find intervals over which the real and imaginary parts of 4>(t) fall within the 
ranges given in the previous section. 



6. Extensions of Theorem Q] 

In this section we discuss the sharpness of Theorem [TJ extend Theorem [I] to a dual family 
of weakly holomorphic modular forms for To (2), and consider analogous theorems for families 
of modular forms for other genus zero subgroups. 

(2) ti 

We have shown that if f k ^(z) is a basis element for the space M»(2) and n is large enough 

compared to \£\, then the majority of zeros of f^(z) lie on the lower boundary of the 
fundamental domain for To (2). We note that the bounds on n are not sharp. For example, 
if £ > 0, then the Fz(u + |) term in the numerator of the integrand takes on its largest 
values near u = 0, while the S^u + |) term in the denominator takes on its smallest values 
near u — ±|, and we have simply taken absolute upper bounds for these functions on the 
interval u G [— |, |] without accounting for interaction between these terms. It is clear that 
sharper bounds on n are possible, and it is natural to ask if the zeros of f^ n {z) always lie on 
z = — 2 + ■ In fact, we now exhibit explicit examples for which the zeros are not on this 
arc. 

Note that if k' = and the degree n + £ = n+ | of the Faber polynomial F(x) is equal to 
1, then we can directly compute that F(x) = x — (8£ — 24). For the root of this polynomial 
to be in [—64, 0], we must have —5 < £ < 3. Thus, if £ > 4 or i < —6, the single nontrivial 
zero of the modular form fk-i+\{z) will not be on the lower boundary of the fundamental 
domain. If k' = 2, the analogous polynomial is x — 8£, and we must have —8 < £ < for the 
zero to be on the lower boundary. Similarly, if k' — and the degree n + £ of the polynomial 
F(x) is equal to 2, we obtain the polynomial x 2 + (48 - 8£)x + (32£ 2 - 188£ + 24), which has 
complex roots when £ < —6 or £ > 6. Such calculations can also be done for polynomials of 
other small degrees to find ranges of £ and n which will not work. Thus, the theorem is not 
true in general, and some condition on the size of n compared to £ is necessary. 

Next, we consider a straightforward extension of Theorem 111 to obtain results for another 

(2) — 

family of modular forms for To (2). Let g k n (z), for n > — I + 1, be the unique modular form 
in M|(2) which vanishes at the cusp and has Fourier expansion beginning 
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These g^ n (z) form a basis for the subspace of M\{2) consisting of forms that vanish at 0. It is 

straightforward to construct g kn (z) by setting g^_ £+1 (z) = S l A F k i(j){z) and obtaining the g£ n 
for larger values of n inductively by multiplying earlier basis elements by powers of <j>{z) and 
subtracting the appropriate terms. These g k n were studied in the k = case by Ahlgren [I], 
and generating functions for arbitrary even weight k were found by El-Guindy jH]- We note 

(2) 

their similarity to the f k n(z). In fact, 

and replacing k with 2 — k and switching r and z shows that the mth Fourier coefficient of 

(2) (2) 

fkn( z ) * s ec L ua l to the negative of the nth Fourier coefficient of g^-kmi 2 )- 

With the above generating function, we can repeat the calculations in Section [3j replacing 
fkn w ith g kn , to see that the quantity e^e~ 7rnsm6 '5'^(— \ + \e ie ) is real valued, and that 
its difference from (— l) n 2cos — vrncos^) is given by an integral almost identical to the 
integral for the fl^(z), except that the </>(r) in the numerator of the integrand becomes a 



4>(z); hence our Hauptmodul expression is the D(z } t) function defined in (4.1). Following 
the argument as before and using the bound for J \D(z, r)| du in Section [4 we find that if 
£ > 0, the difference is at most (.73041) n (71.425)^(l2.50005)(l. 20992), while if £ < 0, it is at 
most (.73041) n (81.38034)^(12.50005)(1.20992). We obtain the following theorem. 

Theorem 2. Let gf^{z) be as above. If £ >0andn> U£ + 7, or if £ < and n > l5\£\+7, 

then at least |_§ + ^"^J of the n + |_f J — 1 nontrivial zeros of g^ n (z) in the fundamental 
domain for To (2) lie on the lower boundary of the fundamental domain. 

Finally, we consider extensions of this method to other subgroups. Consider T (N) = 
{( ™ d) e SL 2 (Z) : c = (mod N)}. If N = 2, 3, 5, 7, 13, then this is a genus zero subgroup 
with two inequivalent cusps and with Hauptmodul 

24 



It is again possible to define analogous families fl n (z) and g\, n '(z) for these groups, and to 
use [8] to find a generating function for each family that mirrors (3.1). Moreover, given a 
modular form / of level N with real coefficients, we find that the quantity e~ l ~ f{— + j^e 16 ) 
is real for < 9 < n; this is a natural place to look for zeros. 

Our argument for the case N = 2 relied on the ability to find a particular horizontal contour 
meeting certain criteria. If we can find a contour passing below the points z = — ^ + -^e ie 
and N * +1 and above all other images of z under the action of T (N), then an application of 
Cauchy's residue theorem to the generating function shows that 

ike _airnsin9 An) f 1 1 %o \ ~ / k9 2nn 2im 



e^e-^W \-- + -e»J - 2co 8 ^ + _ - _ oosO 

_ ikO 27rn sin . - , iri * p * i i * tp 

is equal toe 2 e n times the integral of the generating function along this contour. It 
we take the absolute value of the latter expression, then choosing a contour with horizontal 
height A' means we can pull e~ wn ^N sm0-2A r ) ou tside of the integral, and as long as 4 sin 9 > 
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A', this term has exponential decay as n grows. This means that any bound on the size of 
the remaining integral term proves a result for the location of some of the zeros of fjf^(z) 
for large enough n. 

It is possible to find such a contour for the case N = 3. We take a fundamental domain in 
the upper half plane defined by {z : — \ < $t(z) < \} fl {z : \z — || > |} fl {z : \z + || > |}. 

2-i 



1.5" 




We let k = 6£ + k' for k' e {0,2,4}, and for n > -2£ - [f J we define a basis for Mj(3) 
by letting fj^l(z) be the unique element of M|(3) with Fourier expansion beginning g~ n + 

0(g 2 ^ + ^-l +1 ). We let z = — | + |e* 9 , and require the height of the contour to be above | 
to avoid any images of z other than g^j. We restrict the range of 9 values so that z and 
lie above this contour. The contour may be as close to | as we like, and we may allow 
z to come as close to the contour as we like; once we fix these choices, the integral will have 
a finite upper bound, and the exponential term will eventually dominate when n > Ci for 
some constant C, so that the difference between the weighted modular form fj t n{ z ) an< ^ ^ ne 
cosine function is less than 2. The number of zeros of f^ni 2 ) m ^ ne fundamental domain 
is n + [|J , and as the contour approaches a height of | and z is allowed to approach the 
contour, the number of zeros that can be proved to be on the lower boundary approaches 
.9618n + .2792/c. Thus, we have the following theorem. 

Theorem 3. Let fjfl(z) be as above. If n is large enough compared to £, then the majority 

of the zeros of f^(z) in the fundamental domain for T (3) lie on the lower boundary of the 
fundamental domain. 

A similar theorem clearly holds for the modular forms g^l(z). 

When N = 5,7, or 13, the shape of the fundamental domain is more complicated, and 
choosing an appropriate contour becomes correspondingly more difficult. We leave this as 
an open problem. 
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